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PREFACE 


The tumbling and surface roughness effects on the trajectory 
of entry tektlte are studied in both free molecular and continuum 
flows. It was concluded that, while surface roughness has neg- 
ligible effect on trajectory, the tumbling may play an important 
role In tektlte trajectory and the consequent ablation, provided 
the body shape is different from a sphere . A shape factor B was 
proposed to measure the shape irregularity and was found to be 
a good parameter for correlations between body shape and tumbling 
effects , 
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area 

deceleration 
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Nomenclature (Cont'd) 


S total surface area 

T temperature or total time 

t time 

V velocity 

Vp total volume 

W weight 

angle of attack 
j3 ballistic coefficient 

Y flight path angle 

A shock stand-off distance 

5 boundary layer thickness 

£ density ratio across the shock 

^ speed ratio 

© angle between V and W (Pig. l) 

0 azimuth angle of cylinder axlsj or angle between 

surface normal and velocity 

^ mean free path 

viscosity 

/ density 

to angular velocity 
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Subscripts 




A axial 

C continuum 

e entry condltlonj or edge of boundary layer 

ed end-on position 

FM free molecular flow 

/L roughness 

3 stagnation, or side wall 

sd side -on position 

sp equivalent sphere 

sm smooth 

T transverse 

w wall 

z conditions behind shock 

00 free stream 
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INTRODUCTION 


M 2 ^ 

Most Works of tektlte ablation studies' ^ ' were under the 

assumption that the entry tektlte is of spherical shape and smooth 
surface. The numerical results of those calculations Indicated 
that ablation takes placf. mainly by melt flow rather than by vap- 
orization. This la disturbing in comparison with the observa- 

tions on the recovered tektlte, since melt flow has not been definit- 
ely found on tektlte except for autstrlite and a few Javanltes. 

The effects of surface roughness on boundary layer transition, 
heating and ablation have been studied in Reference 3, which con- 
cluded that the surface roughness on most recovered tektlte was 
formed before earth entry. Consequently, the roughness may trip 
the boundary layers to turbulence and enhance vaporization. Num- 
erical results in Reference 3 Indicated that the effect of surface 
roughness on ablation appeared to be an important factor to ex- 
plain the discrepancy of vaporization vs melt flow. The effect 
of surface roughness on drag has not been discussed in Reference 3 . 

Most of the entry tektltes are expected to be nonspherlcal in 
shape and to carry a not negligible rotational momentum. In 
addition, irregular or asymmetrical bodies, in which the resultant 
of aerodynamic pressure doesn't pass through the center of mass, 
will experience the action of a variable couple influencing rota- 
tion. Thus, one way or the other, the entry tektlte is expected 
to be tumbling in some fashion. Since the drag coefficient is 
highly dependent on the angle of attack, the effective ballistic 
coefficient (or the trajectory) and the total ablation may be sub- 
stantially different from the value of corres poinding spherical 
bodies. In this paper, analysis were conducted under the assump- 
tion that the frequency of tumble is large compared to the change 
in altitude during any cycle; thus, the average drag coefficient may 
be associated with the fixed environment of a single altitude. 

The correction of drag coefficient will change the trajectory of 
entry tektlte, affecting the corresponding albation behavior. 

During earth entry, tektlte experienced various regions of 
atmosphere, from free molecular flow, transition flow to continuum 
flow, depending on the corresponding density level. The division 
of gas dynamics into various regions are usually conducted based 
on characteristic reanges of values of an appropriate Knudson 
number, defined by a dimensionless ratio X /L> where \ denotes 
molecular mean free path and L is a length scale. There are 
various ways in defining X and Lj e.g, the mean free path 
may be measured at freestream or behind the shock, while the length 
scale L may be referred to as the nose radius Rj^, boundary layer 

thickness fl , or the shock stand off distance A, depending on 
the parameters and the environments . 
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Plow In the various transition regimes between free -molecular 
and continuum Is extremely complicated, and no satisfactory 
theoretical solution has yet been offered. No detailed develop- 
ment of drag coefficients or skin friction will be presented In 
this report for transitional flows. Discussions will be concentra- 
ted on the limits of free -molecular and continuum flows. Engineer- 
ing methods will be presented for extrapolating drag coefficients 
and skin frictions between these two well defined regimes. 

In free molecular region, the mean free path Is everywhere 
very much greater than a characteristic body dimension, and the 
distortion of the free stream velocity distribution due to the 
presence of the body Is negligible. However, the actual transfer 
processes are tied up with the fact that the molecules are remit- 
ted or reflected In some manner; thus, a thermal accommodation 
coefficient Is introduced to measure the extent to which the mean 
energy of the remitted molecules Is accommodated toward an energy 
corresponding to the temperature of the wall. The thermal accom- 
modation coefficient is defined by the equation 


a 


e 


Ei - Ea 


( 1 ) 


where Ei. Is the Incident energy per unit surface area per second, 

Ea Is the reflected or re -emitted energy carried away by the 
molecules as they leave the body, and It the energy that 
the re-emltted or reflected stream would have if all the incident 
molecules were re-emltted with a Maxwellian velocity distributions 
corresponding, to the surface temperature T„, Tests of Wledmann 
and TrumplerV^) Indicate that under static”condltlons, the thermal 
accommodation coefficient for air on various typical engineering 
surfaces lies between 0,8? and 0.97. Thus, it is reasonable to use 
unity accommodation coefficient for most practical problems . 

The other important parameter in rarefied gases Is the so 
called speed ratio, t, , which is defined as the ratio of the mass 
velocity to the most probable random speed In the freestream. 

It was estimated in Reference ( 5 ) that the average of air molecules 
at high altitudes to be about 4000-5000 ft ./sec.; an order of magni- 
tude less than the typical speed of tektlte entry. Under the assump- 
tion of unity accommodation coefficient that the air molecules are 
reflected at roughly the same temperature as the tektlte, the re- 
flected particle should have velocities at the neighborhood of 
3000 to 4000 ft. /sec. Since the molecular velocity before and 
after the collision with tektlte Is much less than tektlte velocity, 
the molecular velocity Incident to tektlte can be approximated 
by the tektlte velocity itself and the rebounclng velocity of air 
molecules can be approximated by zero. Thus, the approximation 
of infinite speed ratio, ^ , will be used throughout the analysis 
of free molecular flow . 
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The other limit of the flow regimes Is the well-studied con- 
tinuum flows at sufficiently high Reynolds numbers. The viscous 
effects in this region may be taken into account using classical 
boundary layer theory, with suitable modification for pressure 
interaction if appropriate. Pressure drags on bodies traveling 
with hypersonic speed can be estimated with Newtonian approximation, 
while surface roughness effects can be calculated from classic 
boundary layer theory, With the free molecular and continuum flow 
results calculated, the quantities in transitional flows can be 
interpolated between these two well defined regions. 

Tektite analysis contain a source of error depending on the 
various shape of the tektite, except when the latter is nearly 
spherical such as australite. However, australlte composed only 
a small portion of the recovered tektite, and Indeed most of the 
tektite are Irregular in shape and do not lend themselves to purely 
analytical analysis. Certain plausible assumptions can be made in 
this respect for order of magnitude calculations. To measure the 
geometry of an entering tektite, a shape factor b( 6) is defined 
as the actual surface area divided by the surface area of a sphere 
of equivalent volume. 


B = S/fT (2) 

where S is the tektite surface area and D is the diameter of the 
equivalent sphere, 

D = (6 (3) 

with Vp denoting the volume of tektite. 

Parametric studies for tektite shape effects on trajectory 
and ablation will be carried out In terms of shape factor B with 
circular cylinder used as basic shape for numerical results. It 
is hoped that simple analysis can bring out the order of magnitude 
effect of body shape on tektite entry behavior. 

The surface area and volume of a circular cylinder of diameter 
d and length i are respectively as follows: 


S = f/2 d® + ir dj 

-fh at ('*) 

Substitution of Eq. (4) into (2) and (3) yields the diameter of 
equivalent sphere D and the shape factor B, 

D = 1.1447 d^/^ i 

B = 0.3316 (d/i)^/3 + 0.7631 (i/d)^/3 
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Dependence of B on d/j2 la plotted In Figure 1. It is Indica- 
ted that B tends to Increase in both long cylinder (i/d » 1) and 
short cylinder (£/d < 1) with a minimum at J » d. Both long cylinder 
and short cylinder results will be presented In this report, and 
the inference of these results to tektlte study will depend on the 
actual shape of tektlte. 


TUMBLING AT FREE MOLECULAR FLOWS 


The trajectory of a body entering the atmoaphere from high 
altitude depends greatly on the entry velocity Ve, the entry flight 
path angle 3^^ and the ballistic coefficient 0 which Is defined as 

^-W/OjjA (6) 

Where W Is the weight of the body, and and A are the drag coef- 
ficient and the corresponding reference area respectively. While 
Ve and y'g specify the Initial conditions of earth entry, the bal- 
listic coefficient Is a function of Mach number, Reynolds number, 
and geometry of the entry object. Thus, the key factor for trajector: 
calculation Is the determination of 0* 

Under the approximations of unity accommodation coefficient A^ 

and Infinite speed ratio ^ , the drag on tektlte at free molecular 
flow Is actually equal to the change In momentum of all the air 
molecules striking the surface per unit tlmej l.e. 

- Poc ''J ^ (7) 

Where and A denote the corresponding atmosphere density, 

tektlte velocity and projected front area respectively. Thus, Cj. 
for free molecular flow can be well approximated by 2, based on 
the projected front area. This leads to the following ex press lonr’ 
of ballistic coefficient 


for cylinder at slde-on poaltltlon (cylinder ax5 L flow) 

^ed “ 2W/TTd^ (9) 

for cylinder at end-on position (cylinder axis II flow) and 

^gp = W/(2.0583 d^/3 (10) 

for equivalent sphere, where W denotes weight, and I and d are 
the length and the diameter of the cylinder respectively. 

The same argument leads to the corresponding ballistic coef- 
ficient 


jS* W/(2d£ sin Cl + IT'S d^cosa) (11) 

For a tumbling body, the average ballistic coefficient Is defined 
as being w/Cj^A of an object which has the same deceleration 
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Jh 


aB the average decelerations: 


a 



a(t) dt 


( 12 ) 


where a and t denote deceleration and time respectively, and 

is the dynamic pressure which is assumed not vary greatly during T. 
By definition of average W/C^A 


and, therefore, 




i+i' 




(13) 


( 14 ) 


If the motlon^o_^a cylinder is assume^ tumbling In the trajectory 
plane (l.e. w v = 0 ) and the variation of angular velocity u) 

Is small during T, thus. 


= wt 


(15) 


and the average ballistic coefficient In terms of ot is given by, 



(16) 


substitution of Eq. (11) Into eq, (l6) yields 

5- W/(VfT;^cl + d^) (IT) 

However, since tektltes are probably put out at random spin 
orientation, it will In general be true that 

w • ^4^0 (l8) 

so that the actual will be different from that given by eq. (17). 

Figure 2 shows a {Ji ,"4 , coordinate system with velocity 
vector In the Ic direction, _^e ^vector at an angle 0 with the 
velocity vector and In the k - j plane. The angle^ Is the angle 
between the cylinder axis and the velocity vector v, and ^the 
angle 0 = tJt Is the azimuth angle of the cylinder axls,X^. Prom 
Figure 2 It Is obvious that 


( 20 ) 



cosot mlt * Bln 0 sin 0 ' (21) 

Substituting (21) Into (11)» and (l6) we obta the approxlm:.’te 
solution 

^ (0) » W/[^(2 - i sln^0 - 3/32 sln^0) £d + slnoj (22) 

where a slight adjustment has been made to some of the sine terms 
to allow6(7T/2) matching eq. (17). Thus, at each value of 0 = 0 • 
(0^ *V2) an effect ballistic coefficient >6(0^) can be ° 

determined which represents the average taken over a cycle of 
tumbling In the plane defined by 0 . As Indicated earlier, no 
attempt Is made to predict specific orientation of the plane of 
tumble In the present paper. Instead, an average value of the 
effective ballistic coefficient Is defined for a random orienta- 
tion of the plane of tumbling by the following equation. 

|s(E Ae (23) 

where^ is the probability density. Under the assumption that 
is randomly distributed in space. It can be easily shown that the 
distribution function for © is given by 

P (e<e^) = 1 - cos 0^ 

and the probability density is given by, 

^ - sine 

Substitution of Eq. (25) and (22) Into (23) yields, 

3 =W/( 1.6l66 Jld + 0.7854 d^) (26) 

Equations (8), (9), (17) and (26) represent the ballistic 

coefficients of a cylinder shape at various position and tumbling 
patterns. Their relative values with respect to that of an 
equivalent sphere are summarized in Table 1. Combination of 
Table 1 and eq, (5) yields the numerical results plotted In 
Figure 3, which shows the ratio of ballistic coefficient to that 
of an equivalent sphere as a function of the shape factor B. 

Since the shape of long cylinder differs very much from that 
of the corresponding short cylinder of the same shape factor B, 
the relative ballistic coefficient of an Irregularly shaped tektite 
is expected to be somewhere between or close to those of the long and 
short cylinders . Both results of tumbling In pitch plane and of 
random tumbling are Included In Figure 3 to cover the range of 
various styles of tumbling. Thus, once the volume and the surface 


(24) 

(25) 
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area of a tektlte shape are measured, the range of ballistic 
coefficient can be estimated from the cross-hatched area in 
Figure 3. Since the shape of tektlte can have a shape factor B 
of 2—3» the drag acting on It can easily be 2~3 times larger 
than that of the equivalent sphere. 


TUMBLING IN CONTINUUM AND TRANSITIONAL PLOWS 


As the tektlte approaching continuum flow region, the mean 
free path becomes negligibly small compared to the dimension of 
tektlte and the viscous effects are limited to a thin layer over 
the body surface (boundary layer). Because of the hypersonic 
speed, a modified Newtonian pressure distribution appears to be 
appropriate for mose practical problems. The local pressure is 
given by 

^ “ '’oo ('’s - '’oo> ^ (27) 

Where P Is pressure, 0 Is the angle between the surface normal and 
the velocity and subscripts oo and s denote the conditions at free- 
stream and at stagnation point respectively. The continuity and 
momentum equations across a normal shock are 

^ 00 ^oo ^ ^2 ^2 




Pg + 


V 


Where pand V are density and velocity respectively, and subsc. Ipt 2 
denotes condition immediately after the normal shock. Since 
incompressible flow relations closely approximate actual conditions 
In the stagnation region behind the shock wave, the stagnation 
pressure can be given by 


P 


s 


= Pg + i 




Thus, eq. (27) becomes 
2 = + "oo (2 - 

For Mach number greater than 10, Pqq« and Pg<rs'*^oo “ ^^^ 2 ^ 

and the following simplified relations have been frequently v.sed 
for practical purpose; 

p~<ioo (2 -e) 00^2 

^ Pg COS^ 0 


where g denotes ^ / /*2 * 

The drag coefficient Is given by 


( 29 ) 
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which Is Integrated over the surface of the body exposed to the 
freestream. For the case of a cylinder In side -on position, 

Eq. (29) yields 

= 2/3 (2 - € ) (30) 

which lead the side-on ballistic coefficient 

flsd “ W/ |2/3 (a -6 )| (31) 

in addition, Eq, (29) leads to the drag coefficient and ballistic 
coefficient of the equivalent sphere as follows: 

Oj, = 1/2 (2 -e ) 

^'sp = 1^0.5146 (2 -€) d^/^ 

As to the end-on position, there would be a uniform pressure P 
over the face of the cylinder according to Newtonian theory. ® 
However, Stoney and Swansori(T) found that the average pressure on 
the face of a cylinder 0,909 times stagnation pressure. Using 
this value in Eq. (29) and integrating, we have, for a cylinder 
in end-on position, 


Cj, = 0.909 (a - c ) 

fled = (2 - «) 


(33) 


For a cylinder at angle of attack, the total drag force can be 
approximated by the vector sum of the axial and transverse com- 
ponents 


sin 0t+ F^ cos d 


or 

S* ’co V " % dr °DA (34) 

but 

dr = */oo V = <»oo (35) 

and 

fa> ^A^ = da, oosk (36) 

Combination of eqs. (34), (35) # (36), and the definition of bal- 
listic coefficient yields, 

jg (oL) = W/ ["(2/3 slnotid + 0.714 cosoc. d) (2 -£)1 (37) 
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The average ballistic coefficient for a cylinder tumbling in the 
trajectory plane is given by substituting Eq. ( 38 ) into (I 6 ), 

= W/ j(0.2829 id + 0.3030 d^) (2 - € ) ( 38 ) 

For cylinder tumbling in plane at an angle ( '^/2 - e) with the 
trajectory plane, the average ballistic coefficient can be 
calculated by the substitution of Eqs. (21) and (38) into ( 16 ). 

An approximate solution was obtained, 

fi (e) - W/ [(0.6667-0/48 sin2 e + 0.1 sin^ 8 + 0.303 | sin^ ©) (2-C) <£dj (-J7) 

where a slight adjustment has been made to some of the sine terms 
to allow>5( T/2) matching Eq. (38). The average ballistic coef- 
ficient of a cylinder tumbling in random is given by substitution 
of Eq. (3?) into (23), 

/ = W/ [(0.4 id + 0.1785 d^) (2 -e )] (40) 

Table 1 summarizes the relative ballistic coefficients of a 
cylinder traveling In continuum flow at various positions and 
tumbling patterns. Combination of Table 1 and Eq, (5) yields the 
variation of those relative ballistic coefficients against shape 
factor B as plotted in Figure 4, 

As for drag coefficient in transition reglonspbetween free 
molecular and continuum flows, Bloxson and Rhodes\°' conducted 
a series of drag coefficient experiments in a hypersonic wind 
tunnel on bodies of various shapes with Knudsen number range 
0.0001 - 0 . 34 , The Knudsen number there was defined as I /3 mean 
free path behind shock divided by the shock stand-off distance, l.e. 

Kyj=A/3A (41) 

The drag coefficients so obtained on spheres are plotted 
in Figure 5 as a function of Knudsen number, which shows that drag 
coefficients can be presented as an exponential function of Knudsen 
number in transition flow. 

The most remarkable point about Figure 5 is that all of the 
shapes considered have the same variation with in the transi- 
tion region and the same dra^T coefficient in the free molecular 
region. The variation seems to hold for any shape and gas, having 
only the requirement that the flow is over Mach number 4. An 
empirical formula was derived b£ sed on the data of spherical 
bodies in Reference 8; i.e. 


1.52 



where C^, and Cjj, ^ denote the drag coefficients In free- 
molecular Id continuum regions respectively. 
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ROUGHNESS EFFECT 


Surface roughness Is known to promote laminar-turbulent transi- 
tion and to Increase surface drag and heating by a significant 
amount . The surface drag on rough surfaces depends on the number 
of roughness elements per unit area, on their shapes and heights, 
and on the way In which they are distributed over the surface; 
since the main contribution of surface drag la from the pressure 
force acting on the roughness elements. The number of parameters 
required for roughness description Is extraordinarily large and no 
analysis have been developed to include all the parameters required. 
Most authors measured the roughness effect on surface drag by a 
single parameter of effective roughness defined as the rough- 

ness size of closely packed sand roughness which gives the same 
value of skin friction coefficient as the actual roughness. 

In this way, the difficulty has been shifted to the deter- 
mination of effective roughness from the measurement of actual 
roughness. However, for order of magnitude analysis such as 
tBktlte studies, the average height of the reasonably spaced rough- 
ness elements can be used as effective roughness for practical 
calculations. As discussed In Reference 3, an effective roughness, 
Keffi equal to 40 mils Is assumed for the tektlte studies in this 

report . 

The exact way to assess the surface roughness effect on 
trajectory Is to calculate the drag coefficient of a roughened 
tektlte which Is tumbling during entry. However, the coupling 
phenomena of roughness and tumbling is a very complicated problem 
because of the cross flows Involved In an Inclined cylinder, which 
is beyond the scope of this report. The approach here Is to 
estimate the roughness effect according to the most noticeable 
case, l.e,, the case that roughness has maximum effect. For a 
long cylinder (d/i<l), the roughness effect Is maximum as the 
cylinder Is at end-on position, and the numerical results of 
roughness effect will be carried out In this position. 

Roughness effects on drag are closely related to the behavior 
of skin friction drag coefficient through various regions of 
earth entry. Reference 9 presented an empirical relation of the skin 
friction coefficients in the form of (C~ M) as a function of the 
interaction parameter III/ VRje, which Is ^defined for flows with 
significant viscous effect. In such a case, the significant 
characteristic dimension of the flow field Is the boundary layer 
thickness 5 rather than a dimension L typical of the body itself, 
and the corresponding Knudson number based on the length scale 
Is actually proportional to the Interaction parameter M/ ^e. 
Hoerner(9) suggested that for skin friction coefficient simple 
gas dynamics for continuum flows are applicable to m/ i/^< 10 “ 2 , 
and a transitional phase between 10“^ and 10+3., above which free 
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molecular flow 1b finally established. Figure 7 presented both 
laminar and turbulent skin friction coefficient in terms of 

vs m/ ^Re. Very good correlation is found between free molecular 
and continuum results , Naturally, turbulent friction should not 
be expected to reduce belcw the level as found for laminar bound- 
ary layer flows. The experimental results as plotted in Figure 6 
suggout that the laminar function represents the limiting condition 
for the turbulent function. It appears that a low Reynolds 
number the distinction between Igimlnar and turbulent boundary 
layer flow becomes lost. It is suggested^} that the vorticity 
originating across a boundary layer induced shock wave behaves like 
turbulent, thus possibly rendering a ’’laminar" layer effectively 
turbulent. For surfaces containing roughness of 40 mils such as 
observed on many recovered tektite, the boundary layers will be 
tripped to turbulence somewhere between m/ ss 10"1 and 1 
according to the numerical results in Reference 3; and turbulence 
exists in boundary layers in most of the transitional region. 

Thus, it is plausible to use the conventional analysis of rougn- 
ness effect in turbulent flow through the whole trajectory of 
tektite entry. 

The effect of roughness on turbulent skin friction has re- 
ceived considerable Investigation in low speed flow as summarized 
by Clauser,^^®/ The roughness skin friction in high speed flow was 
studied by Goddard, who showed that the effect of compresslblillty 
can count as a reduction in wall density as Mach number increased. 

In addition, the skin friction Increase due to roughness correlated 
with the roughness Reynolds number Rej^*s ^ Ur K/M w 

Independent of Mach number. The data of Reference 11, which were 
limited to adiabatic flow conditions, yield a close approximation 
for roughness effect on skin friction: 

( 12 ) 

More recently Young' ' found that as the wall temperature was 
reduced, the skin friction dropped below the Goddard value for 
adiabatic walls. Reasoning that the effect of heat transfer was 
to produce a density gradient near the wall such that the effective 
density at the roughness surface was less than the wall value. 

Young developed a reference temperature method that correlated his 
results, 

Cf/Cf^= 0.365 (Te/Tr) f 0.635 (Te/Pw) (44*) 

where is the incompressible value of including roughness 

effect ^ 

The relation between roughness skin friction coefficient and 
its adiabatic wall value can be derived directly from Eq. (44) 

i.e. 
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0.365 (TwAr) + 0.635 


(45) 


Combination of Eqs, ( 45 ) and ( 43 ) yields the relation be- 
tween rough wall skin friction Including heat transfer to that of 
smmooth wall under adiabatic conditions \ ^3) 

Cf/Cf = |o .365 (TwAa) + 0.635] log 10 Rej^* (46) 

Thus, once the reynolds number around the roughness elements, 
the roughness size K and the temperature ratio Tw/Tr are given, 
the affect of surface roughness on skin friction can be calculated. 
The skin friction coefficient of smooth wall can be estimated 

from the correlations In Figure 6 . ° 

The total drag for a long cylinder traveling In an end -on 
position can be expressed as follows; 

= '*co ^ ^ + C ) (4-7) 

where term denotes the friction drag acting on the sidewall of 

the cyllfider as Indicated In Figure 6 and Cp term Indicates the 

pressure drag acting on the end face of the cylinder. Usually the 
friction drag for a smooth cylinder Is at least three orders of 
mamltude less than the pressure drag, providing the area rateo 
Ag/Ag Is of order one . This explains why only pressure drags are 

considered In tumbling studies. For rough surfaces, the friction 
drag can be larger than that of smooth wall, and Its effect can 
be expressed by the relative drag of rough to smooth long cylinders 
as follows; 



r-sm 


F 


D,r 

D,sm 



(48) 


Where Cp can be estimated from Eqs. ( 9 ), (33), and (42) for free 
molecular, continuum, and transitional flows respectively, 

^0 

from Figure 6 , from Eq. (46) j from hypersonic results 

such as in Reference 3, and A^/A^ from geometry considerations. 

The area ratio A_/a. can be related to the shape factor B through 

s s 
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Eq. (5). 

Figure 7 shows the relative drag of rough long cylinder to 
smooth one at end-on position, In which drag Is presented as a 
function of shape factor B. The Mach number here Is a constant 
value of 30, which Is the average value In the actively ablating 
region of tektlte entry, (3) 

The Reynolds number here, based on the equivalent diameter 
of one inch, is 105, which is the upper limit that an entry 
tektlte may experience according to the numerical results In 
Reference 3. The ratio of -roughness size to the equivalent diameter, 
I^D, is selected to be 0.04j since the average roughness on tektlte 
is about 40 mils and the average diameter of tektlte is probably 
at the order of one inch. Figure 9 indicates that the drag In- 
crease due to roughness at Re^ = 105 Is less than 10^ for usual 

tektlte geometry (say KBoT* Since roughness effect on drag 
Increases with Reynolds number, the drag Increase indicated In 
Figure 7 gives the upper limit of roughness effect on drag for 
tektlte entry. The Reynolds number of an entry tektlte usually Is 
much less than 105; thus. It Is plausible to assume that the 
surface roughness effect on drag can be neglected for tektlte 
entry studies . 


TRAJECTORY AND ABLATION 


The trajectory of a vehicle entering the atmosphere from great 
altitude depends greatly on the entry velocity the flight path 

angle / , and the ballistic coefficient's , In general^ both the 
flight path angle and the ballistic coefficient vary due to the 
presence of gravity force and to the change of drag coefficient 
during flight. Hence, It Is clear that exact solution of the traject- 
ory Is formidable. However, In the trajectory region of active 
eolation, the velocity of the entry tektite is so high that its 
total drag is much larger than the gravity force and that the drag 
coefficient Is nearly constant because of the corresponding high 
Mach number environment . v / Thus, for a simplified trajectory 
analysis like this. It is appropriate to assume that both K* andj9 
are constant. In addition, the atmosphere density/® , which 
is a function of altitude only, is given by the rel^ion,. 

where and X are constant equal to O.O9I lbs. /ft. ^ and 23120 ft, 
respectively according to the curve fit for high altitude atmosphere 
In Figure. 8. Based on these assumptions the relations, Allen and 
EggersV^^^ derived the following relation 

V = Vg exp r./>^ A/ (2 )j (50) 

Which can be used to estimate the vehicle velocity as a function 
of altitude once the parameters of the entry conditions, such as 

j3 ^ Vg, are given. 

Numerical results In Reference 3 Indicate that ablation rate 
was practically zero if velocity reduced to or below 10,000 ft/sec. 
Thus, one can define a consolidating altitude at which the velocity 
Just reaches 10,000 ft/sec and no significant ablation there- 
after. Combination of Eqs. (45), {SO), and (5) yields the 
numerical results plotted in Figure 9 which show the consolidating 
altitude as a function of shape factor B, The ballistic coefficient 
formula used In Figure 9 is based on that of random tumbling In 
free molecular flow with the equivalent diameter of one Inch, 

It Is clearly indicated that the shape factor B has significant 
effect of tumbling on trajectory. In addition, the geometry style 
(long or short cylinder) has only slight effect on the consolidating 
altitude. Indicating that shape factor B Is a fairly good para- 
meter to describe Irregular shapes. Figure 10 shows the similar 
conclusion for the trajectory calculations based on the randomly 
tumbling ballistic coefficient In continuum flows. 


For an order of magnitude estimation, the total recession of 
an entry body can be measured by the total heat transfer to the 
body. The stagnation point heat transfer ra^-' to a sphere of 
be approxlmlmately calculated the following 


radius R«.caD 

equation” (15) 


where the unit of heat transfer rate q„ is In BTU/ft^-sec, radius 

a 

In feet, and velocity V in ft/sec. Since the most obvious 

effect of tumbling on tektite behavior Is velocity reduction due 
to higher drag. It is plausible to measure the ablation rate by 
the substitution of reduced velocity Into Eq. (5r/). The relative 
total ablation of velocity reduced tektlte with respect to that 
of equivalent sphere can be measured by the total heat transfer 
to the stagnation point, i.e.. 


A^tumb ^^tumb^ 




whereA S Is total recession, subscript tumb and o denote tumbling 
correction and no tumbling respectively. Combination of Eqs . (ST), 
(Si), and (5) yields the numerical results plotted In Figure 10, 
where the relative total recession Is presented as a function of 
shape factor B. The ballistic coefficient formula used in Figure 11 
Is based on that of random tumbling In free molecular flow with 
the equivalent diameter of one Inch, The shape factor B has 
significant effect on total recession, while both / and geometry 
style (long or short cylinder) have only negllgible®effects . 

Figure 12 shows the same conclusion for the total recession cal- 
culations based on the ballistic coefficient on continuum flow. 
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CONCLUSION 


This study indicates that surface roughness effect on the 
total drag of entry tektlte is negligible, and can be Ignored 
for tektlte studies. For shapes other than sphere, the tumbling 
of tektlte has significant effect on the total drag and on the 
consequent entry trajectory. In addition, the shape factor B 
was found to be a good parameter for the measurement of Irregular 
shape and its tumbling effect on tektlte entry. Rough estimation 
of the tumbling effect on the total recession can be estimated 
once the shape factor B and the size of tektlte are given. 
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ABLE I 


ReUtive Ballitt Coefficient of e Cylinder with Respect 
to the Equivalent Sphere (/^aphere^ /^cylinder ^ 



Flow Reaime 


Caee 

Free Molecular 

Continuum 

Side -on 
End -on 

f'.S'Tl? 

d 2/3 
0. 7632 if 

£. 1/3 
1. 2955 (^ 

d2/3 
1. 3875 (j) 

Tumbling in 
pitch plane 

0.6186 (J) + 0.4858 (2) 

0 1/3 d 2/3 

0. 5497 (i) + 0. 5888 (j) 

Random 

tumbling 

^1/3 d2/3 

0,7852 (J) +0. 38l6(j) 

^1/3 d2/3 

0.7773 + 0. 3469 (p 

d ^ 
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FIGURE I - VARIATION 
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FIGURE 3 - TUMBLING W , IN FREE MOLECULAR FLOWS 



25. 



FIGURE 4 - TUMBLING EFFECT AND SHAPE FACTOR IN CONTINUUM FLOWS 






FIGURE 5 - VARIATION OF DRAG COEFFICIENT IN TRANSITION REGION 




RACTION 
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FIGURE 12 - TUMBLING EFFECT ON TOTAL RECESSION (CONTINUUM) 



